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Abstract. We investigate random, discrete Schrodinger operators which arise 
naturally in the theory of random matrices, and depend parametrically on Dyson's 
Coulomb gas inverse temperature f3. They are similar to the class of "critical" random 
Schrodiner operators with random potentials which diminish as | a; | 2 . We show that 
as a function of (3 they undergo a transition from a regime of (power-law) localized 
eigenstates with a pure point spectrum for [3 < 2 to a regime of extended states with 
singular continuous spectrum for (3 > 2 . 



1. Introduction 

Dyson's Coulomb gas model for the spectral fluctuations of random matrix ensembles 
was recently formulated in terms of ensembles of symmetric, real tridiagonal matrices 
PQ; see also [21 El Ej. These ensembles share the property that the diagonal matrix 
elements are independent, identically distributed Gaussian random variables, while the 
off diagonal elements are independent random variables whose probability distribution 
function (PDF) depends both on the position within the matrix and on the inverse 
temperature j3. We consider these matrix ensembles as ensembles of discrete Schrodinger 
operators, with random on-site potentials (diagonal matrix elements) and random 
hopping amplitudes (off-diagonal elements) with prescribed PDF. This viewpoint has 
been the starting point of recent studies into characterizing the largest eigenvalues of 
the limiting matrix ensemble in terms of a certain stochastic Schrodinger operator [HIE]. 

The interest in these operators stems also from the fact that they are analogous to a 
class of operators for which the random potential diminishes as a power law where 
x marks the position along the chain. Similar systems were thoroughly discussed in the 
mathematical literature (see e.g, [7J [8j [91 10]), where it was proved that a decaying 
diagonal disorder with a < 1/2 induces localization and the spectrum is pure point. 
However, for a > 1/2 the states are extended and the spectrum is absolutely continuous. 
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The behavior at the critical power a = 1/2 depends on the details of the potential, and 
the eigenstates can be either power-law localized or extended. The model we study here 
is related to this critical class, but not exactly, since in the present case the transition 
amplitudes are also random variables. We show that in this model, the parameter (3 
determines the spectral properties: in the regime < (3 < 2 the spectrum is pure - 
point, and the eigenstates are power law localized, while for (3 > 2, the eigenstates are 
extended and the spectrum is singular continuous with a (3 dependent spectral measure 
with a Hausdorff dimension \ — \. 

We start with a short survey of the relevant information from Random Matrix 
Theory (RMT). The random matrix ensembles GOE, GUE and GSE are ensembles 
of N x N real symmetric, complex hermitian or hermitian real quaternion matrices, 
respectively, whose matrix elements are independently distributed random Gaussian 
variables with joint distribution proportional to 

exp(— c Tr H 2 ). 



The probability distribution functions of their eigenvalues Ai, 
a concise form 



(1) 

A 7v can be written in 
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Here, (3 = 1,2,4 is used for the GOE, GUE, GSE ensembles, respectively, Gpx are 
known normalization constants, and c in ([1]) has been chosen to equal 1/2, 1/2, 1/4 for 
the GOE, GUE, GSE, respectively. It can be shown from (j2J), or alternatively directly 
from the definitions of the ensembles by studying their resolvent, that to leading order 
the normalized spectral density is supported in the interval [—\/2{3N , y/2(3N] and it 
assumes the "semi-circle" law 



P(A) 
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A 2 
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Recently, a systematic way to construct the ensembles corresponding to arbitrary 
(positive) (3 was introduced [I], and it is based on the following observation [TTI [12l [13] . 
Any real symmetric matrix A G GOE can be orthogonally transformed to a tridiagonal 
form 
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The probability distribution function of the matrix elements of the corresponding 
tridiagonal matrix TCn has the following properties: 
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• The diagonal elements {a n } are real, independent, identically distributed, Gaussian 
random variables. 

• The off-diagonal elements {b n } are non-negative, independently distributed random 
variables, with PDF 

Pooe(K) = Xn{b n ) = f ^y(^) n " 1 e" fe " • (5) 

The surprising new result is that by distributing the off diagonal matrix elements using 
the PDF 

Pp(K) = X n(b n ) = -irWH""* , (6) 

V 2 / 

the eigenvalue PDF of Hn is given by (J2]) for any positive (3. Thus, the study of 
the tridiagonal ensembles (denoted by G(3E) provides a convenient way to interpolate 
between the classical random matrix ensembles with the discrete (3 — 1, 2, 4. (A similar 
method was recently applied in [2] to Dyson's ensembles of unitary matrices, to get 
Circular (3 Ensembles; this was used in [H] to calculate eigenvalue statistics for CMV 
matrices). 

Denoting by the expectation value with respect to the G(3E measures, we can 
easily find, 

and, 

((bn-(bn)) 2 )p=\ + 0(^j . (8) 

Thus, for large n, the PDF ([6]) limits to the Dirac distribution 5{u n — 1) in the normalized 

variable defined by b n = \J~^u n . This also shows that by scaling the matrix elements 

H N i-> ^/2/PNHn, the new off diagonal elements decay as n ' where n is counted 
from the bottom row of the matrix. 

Once the matrix A under consideration is in tridiagonal form (j3J), a simple recursion 
relation can be written for the characteristic polynomial D^{\) := det (XI — A) = 
det (A/ — Hn)- Denoting the determinant of the top n x n sub-block of XI — Hn by 
D n (X), expansion by the last row shows 

D n = (A - a n ) ZVx - b 2 n _ x £> n _ 2 ; 1 < n < N , (9) 

subject to the initial conditions 

D_i = 0, D = 1 . (10) 

We remark that by computing the zeros of the characteristic polynomial for the 
tridiagonal matrices (jl]) one is sampling from the correlated PDF (T2l) . 

The matrix H.n, in the limit iV — > oo can be considered as a representation of 
a discrete quantum hamiltonian which governs the dynamics of a quantum particle 
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hopping randomly between sites on the half line. The distribution of the "on-site 
potentials" and "hopping amplitudes" are provided by the PDF of the a n and the 
b n respectively. In the mathematics literature, this is referred to as a discrete random 
Schrddinger operator, or a random Jacobi matrix. We address the following questions: 
i. Whether, for almost all realizations, the eigenfunctions of the random hamiltonian 
are localized, or in other words, if the spectrum is continuous or discrete, ii. In what 
way the localization depends on the parameter (3. 

Consider the matrix (TjO) for a finite N. The eigenvectors v = {y\, • • • , Vn) satisfy 

U-w = Av bn^Vn-t + (a n - X)v n + b n v n+1 = , V 1 < n < N , (11) 

with the boundary conditions 

v = v N+ i = . (12) 

The homogenous boundary conditions ( fl2l) can be satisfied only for N discrete values 
of A, and this set coincides with the zeros of the characteristic polynomial Pat(A). 

Ignoring the boundary conditions for a while, the recursion relations ( Till can be 
solved for any A and N, by two independent vectors x and y. The Wronskian 

W(x, y) = b n (x n+1 y n - x n y n+1 ) (13) 

is independent of n and therefore the growth rate of the vectors compensate each other 
so that the Wronskian remains constant. It is convenient to chose one of the vectors, 
say x as the solution which satisfies the initial conditions, 

x = 0, xi = 1. (14) 

It can be computed (for any A) by forward iterations of ( TTTT) . Comparing the two initial 
value problems AS]) . (1T21 and (ITT1), (ITi]). we find that 

n 

D n = x n+1 Y[ b m , (15) 

m=l 

which can be proved by direct substitution. Note that the forward iterations usually pick 
up the solution with the fastest growing rate. An independent solution of the recursion 
relation can be obtained by imposing the condition yjv+i = 0, y^ = 1 at an arbitrary 
value of iV and perform a backward iteration of (TTTl) . This can be done for every A, 
and in most cases the solution to be picked up is the one for which \y n \ is the fastest 
increasing solution when n is decreasing (for < n < N). The Wronskian relation ( JTBl) 
implies 

bj^XN+l 

Vo = —r • (16) 

°o 

Before addressing the effect of randomness it is useful and instructive to study first 
the one parametric family of mean hamiltonians which are obtained by replacing a n 
and b n in (T4]) by their G/3E expectation values. This way we can better appreciate the 
effect of randomness on the quantum dynamics. We shall show that the eigenfunctions 
of the mean hamiltonians are extended, and the spectra are absolutely continuous for 
all p > 0. 
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The mean hamiltonians (TC)p are tridiagonal matrices with vanishing diagonal 
matrix elements. The off diagonal terms are given by ([7j), and, to leading order, are 
proportional to y/n. Thus, for large n, the recursion relations for the components of an 
eigenvector are: 

y/n — lx n -i + Vnx n+1 = V2 Xx n , (17) 

where A = The solution of this recursion relation subject to the initial condition 
Xo = 0, xi — 1 can be written in terms of the normalized eigenf unctions of the one 
dimensional harmonic oscillator 

/ 1 \^ _a£ 

X n+1 = U n (X) = f r^ n[2 n J e ~ H nW , (18) 

with rt_ 1 (A) = 0. The completeness and orthonormality of the Hermite polynomials 
implies that for any real A, /i, 

oo 

J~] u m (X)u m (ii) = 5(X - //). (19) 

m=0 

This proves that the spectrum of the operator (TCn)p for iV — > oo is absolutely 
continuous and supported on the entire real line, for all /3 > 0. For finite matrices, 
the boundary condition Ujy+i = is satisfied if A is chosen as one of the zeros of the 
Hermite polynomial H^(X). For finite but large N the spectrum is located in an interval 
of size 2\/2N centered at A = 0. The normalized spectral density p(/i = X/y/2N) is 
supported on the interval [—1,1], and approaches the semi-circle law 

p(ji) = -xjl-p 2 (20) 

7T 

in the limit iV — > oo. 

In the subsequent paragraphs, we shall show that, in contrast with the 
eigenfunctions of the mean Schrddinger operators which are delocalized, the 
eigenfunctions of the disordered operators are power law localized for the G(3E ensembles 
with (3 < 2. Beyond the critical value (3 = 2 the eigenfunctions of cannot 
be normalized and the spectrum is continuous. However, the disorder has the effect 
that now the spectrum is singular continuous with a spectral measure which has a f3 
dependent Hausdorff dimension 1 — i . 

A prominent quantity of interest in the study of random Schrddinger operators is 
the mean growth rate of the eigenvectors x [T7] . It is related to the properties of the 
characteristic polynomial by 



C p = -( log 

n 



X 



n+1 



1 1 1 n 

(\og\x n+l \) = (log\D n \) p + - VVlogl&J^ . (21) 

n \ Ian m ' * 



n \ I p n n 

m=l 



Thus, the mean Lyapunov exponent Cp which characterizes the Anderson model, is 
expressed in terms of the expectation value of the logarithm of the characteristic 
polynomial of the Gj3E ensemble. Since the latter is known from RMT, and the mean 
value of the rightmost term in (|2T|) can be evaluated directly, the mean Lyapunov 
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exponent for this model can be written down for any value of A. Using the exact PDF 
for the b n , we get 



This is derived starting with the identity [16] 

N-l c N-l 

I] r(a + 1+jc) = c ^(iV-l)/2+iV(a + V2) (27r) -iV(c- 1 )/2 JJ JJ + J 

3=0 p=l 3=0 C 

valid for c G Z+. Let 

f u ^ - ^" TT G{n + {a - p)/c ± V g (-p/ c + ^ 
M ' J Ji G(n-p/c + l)G((a- P )/c+l) ' 

where G denotes the Barnes G-function. Using the asymptotic formula 
(G{N + a + l)\ „ WT . a-6, „ A , WT .a 2 -^ 



log U(iv + 5 + i) J *~ J 6 - + — log 2ff + l (a - 6)iV + — J log N + o(1) 

it follows that 

f n (a, c)~ exp(an log n)c^e— n -(-iW*+*V* JJ ^(-p/c + 1) 

Taking logarithms of both sides, differentiating with respect to a, then setting a = — 1+c 
gives (1221) . Furthermore, if c is rational, c = s/r for s,r£ Z + , the identity [16] 

r-l 



f fa s /r) = Y] fn{a + SU/r,S) 



v=0 



proves that (j221) remains valid for n a multiple of r. Thus if the limit is to exist for 
non-integer c it must be given by (1221) . 

The G(5E expectation value of log |D n (A)| is given by 



i(log|£> B (A)|> /9 = iV<log|A-A 1 |> /J 

n n 

i=i 

= J dypp(y) log \y - A| = ± (log ^ - l) + O {^\ , (23) 

where the exact spectral density was replaced by its semi-circle limit Q, and A <C y/2n/3. 
Substituting in ( 12T1) . we find that 



ifi-il 



logn 2V2 f>> 

= — • 24 

n 

Thus, on average, the components of the eigenvectors x behave as 

|x n | 2 xn^. (25) 
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Using (TIB"]) we expect the other solution of the recursion relation to exhibit a mean decay 
rate of 

|y n | 2 xrT?~i (26) 

This result suggests the following scenario: The eigenvectors are square normalizable 
when (3 < 2 which would imply that the spectrum is pure point in this (3 domain. A 
transition to a continuous spectrum would be expected in the complementary domain. 
Indeed, with a little more effort, we could show that this is true. 
Consider the matrices: 

( A— q n _ K-i \ 
K K \ (27 ) 

and their product 

rpX qA qA qA (OQ\ 

1 n — °n ' °n-l ' ' ' °1 ' 

T£ have the property that for the eigenvectors x 

+ l 1 /; A I I- (29) 



x 



n 




In fact, the above holds for any vector satisfying equation (fTTj) (with any boundary 
conditions). Thus, in order to control the asymptotics of x, it is reasonable to try and 
control the asymptotics of || ||. By an adaptation of methods from [W\ (for details 
see [IS]), it is possible to prove 
Proposition 1. For any A £ 1 

log II T A || 2 1 1 
lim ° g " = - - - (30) 

n^oo logn [3 2 

with probability one. 

With this at hand, a direct application of the methods of [10] (see also [8j) imply 
Theorem 1. For any A G R, with probability one, equation (11) has a solution y 
satisfying 

\y n \ 2 x rr ( 5 + K (31) 

Any solution to equation (llip that is linearly independent from y, satisfies 

|x n | 2 xn^"i (32) 

The asymptotic behavior of eigenfunctions is intimately connected with the spectral 
measure of H, associated with the vector 5i = (1,0,0, . . .). This is the object defined 
by 

N 

d/i(E) = w - limjv^oo ^ I ( s i I VVi) f $(^- K) ( 33 ) 

n=l 

where X n are the eigenvalues of H-n, and ip n are the corresponding normalized 
eigenvectors, using the technique of spectral averaging (see e.g. [20J), the following 
can be shown to ensue from the theory of subordinacy 
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Theorem 2. For any (3 the essential spectrum of H is R. 

If (3 < 2, then, with probability one, \x is pure point with eigenfunctions decaying 

as 

\x n \ 2 x n~ { ^\ (34) 

If (3 > 2, then with probability one, for any e > 0, /i is absolutely continuous 
with respect to (1 — | — e:)-dimensional Hausdorff measure and singular with respect to 
(1 — | + e)-dimensional Hausdorff measure. 

Thus, we see that, as long as (3 < 2, H has square-summable eigenfunctions whereas 
for (3 > 2 the spectrum is purely continuous. This spectral transition at /3 = 2 from 
pure-point to continuous spectrum is, in a certain sense, continuous in /3, since the decay 
rate of the eigenvalues (| + -|) changes continuously in /3. This is in striking contrast 
with the transition expected in the Anderson model. Thus, for the case studied here, 
the Inverse Participation Ratio, for example, should change continuously from 1 (when 
(3 = 0) to (for (3 = oo). 

An interesting feature of this continuous transition is the connection between the 
'extendedness' of states and the level repulsion observed on finite scales. Equation (2) 
above implies that the probability of finding pairs of close eigenvalues, of the finite 
dimensional matrix 7{n, diminishes as (3 increases. Theorem 1 says that as (3 increases, 
the states decay at a slower rate. These two facts are complementary: We expect slower 
decay rate to be connected with stronger level repulsion. Our analysis confirms this 
expectation and gives physical meaning to this aspect of the eigenvalue statistics of 
G(3E. 
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